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Spectrum of SU(2) SUSY Yang-Mills Theory with a light gluino 
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We report on new results for the low lying spectrum of N=l SUSY Yang-Mills Theory with SU(2) as the gauge 
group. Simulating on larger lattices at k = 0.194 and (3 = 2.3, we slowly approach the supersymmetric limit at 

^gluino — 0. 



1. Introduction 

Progress has been made in further determin- 
ing the low lying spectrum of N=l SU(2) SUSY 
Yang-Mills Theory (SYM). The motivation for 
doing large scale simulations has been presented 
at previous conferences and in various publica- 
tions, for a review on the subject see Q. 

This work is related to a previous project of 
the DESY-Miinster-Roma collaboration to simu- 
late SYM in the vicinity of the supersymmetric 
point |l| using the dynamical-fermion two step 
multibosonic algorithm (TSMB). This involves a 
light gluino mass. Previous calculations left open 
the questions a) how close to the SUSY limit was 
the model actually simulated and b) the multi- 
plet structure of the spectrum. The first question 
was answered by last year's SUSY Ward-Identity 
results p| showing that the gluino was heavier 
than expected from previous estimates. The sec- 
ond question is still under investigation and up to 
now not conclusively answered. Here we attempt 
an analysis of the spectrum closer to the SUSY 
limit, i.e. with a lighter gluino, where the SUSY 
pattern of masses should become more apparent. 



2. N=l SYM on the lattice 

A basic assumption about the non-perturbative 
behavior of N=l SYM is confinement, as in QCD. 
Following results from effective action analyses 
(I), we expect to see two chiral multiplets at 
the bottom of the SYM mass spectrum. Their 



content is a spin-i fermion (the gluino-glueball) 
and two bosonic states with opposite parity: a 
scalar and pseudoscalar gluino-gluino bound state 
(gluinoball) a — fo and a — rjf (a denoting the ad- 
joint representation of the gluinos) and the scalar 
and pseudoscalar glueball with J PC = h and 
0++. 

On the lattice, Poincare invariance is broken 
and therefore SUSY. We use Wilson fermions 
where SUSY is also explicitely broken by the lack 
of chiral symmetry. Finally a soft breaking is 
caused by the gluino mass. Because of all this 
SUSY breaking the SYM spectrum gets distorted 
on the lattice in an essentially uncontrolled way. 

We compute the SYM mass spectrum on the 
lattice from first principles using the familiar 
techniques known in QCD. The correlation func- 
tions of the gluinoballs are similar to those of 
QCD flavor singlets 

C S g(At) = (Tr sc [rA x jTr sc [rA w ] 

X 

-2TT sc [TA xy TA yx }) . 

where y is a fixed source, the trace is over spin 
and color and T e (1,75). 

The gluino-glueball is associated in the con- 
tinuum with the interpolating operator = 
c^yTrcfF^yA] where X(x) is the gluino field; on 
the lattice the correlation is 
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where the trace here is only over color. The 
glueball operators are the equivalents to those of 
QCD. 

In order to be able to make meaningful state- 
ments of the ermergence of SUSY in the contin- 
uum through extrapolations, we want the gluino 
to be as light as possible. As in QCD, the natural 
barrier is algorithmic performance due to critical 
slowing down. 

3. Numerics 

Following Curci and Veneziano, we employ 

Seff[U] = 0J2 (l - ^ReTWp)) - I logdetQ 

(Q is the usual fermionic matrix) from which we 
see that the gluino effectively has flavor number 
Nf = i . Expectation values of operators read 

(O) = Z- 1 J V[U]TT[M]0{U)e- s °^ [u] 

where VT\M\ is the Pfaffian of the antisymmetric 
fermion matrix M = CQ. 

The most suitable algorithm for this model is 
the two-step multiboson (TSMB) algorithm 0]; 
for details on its application to SYM see M. It 
relies on representing the fermion determinant in 
the form 

\det(Q)\ Nf - m (2] . ■ 

The polynomial approximations satisfy 

pW(g 2 ) ~ X - N fi* 

lim P£(x)Pif(x) = x~ N f'\ x€[e,X] 

ri2 — >oo 

where the eigenvalues of Q 2 = Q^Q on a typical 
gauge configuration are required to be in the in- 
terval [e, A]. Pni gives a crude estimate of the 
fermionic measure and is used in the bosonic rep- 

(2) 

resentation of the determinant. _P ri2 is a cor- 
rection factor that is accounted for by a global 
accept-reject step. The algorithm is made exact 

(3) 

by a third polynomial P„ 3 through reweighting 
the gauge configurations in the expectation val- 
ues. 



For the results presented we use the following 
samples, all at (3 = 2.3: 



K 


LxT 


€ 


A 


Stat 


0.1925 


12x24 


0.0003 


3.7 


4204 


0.194 


12x24 


0.0001 


4.5 


2034 


0.1955 


12x24 


0.0000125 


5.0 


5324 


0.194 


16x32 


0.0002 


4.0 


664 



where we also indicate the e and A used in the 
simulations. The configuration on the 12 3 x 24 
lattice were produced in || and 0. In order to 
check finite size effects we started a new produc- 
tion run on a 16 3 x 32 lattice; the configurations 
on the last line of the above table are a snapshot 
of it. 

From the SUSY Ward-identities study g] it fol- 
lows k c (/3 = 2.3) « 0.197, where k c (P) is the hop- 
ping parameter where the gluino is massless. 

4. Results 

The new results presented here essentially re- 
fer to the two samples of configurations on the 
12 3 x 24 lattice, k = 0.194, 0.1955 produced in 
These correspond to a lighter gluino com- 
pared to the previous extensive analysis of the 
spectrum accomplished in Q where the largest k 
value was 0.1925. We also present preliminary re- 
sults for the bigger 16 3 x 32 lattice at K = 0.194. 
The available results are summarized in table |l|. 

Glueballs. We used APE smearing with 
N A pe = 36, 40 and e APE = 0.285 for the glueball 
operators. We then averaged the effective masses 
over the emerging plateau. The data for the h - 
glueball is generally very noisy so that averages 
can only be taken at time-slice sparations 1 and 
2. The case of the ++ -glueball is more favorable 
(see fig. |l|). The statistics for the larger lattice is 
still too low to determine a mass. 

Gluinoballs. The connected part of the correla- 
tion function was evaluated by choosing a random 
source, the disconnected by using the volume- 
source technique, both without smearing. The 
signal for the a — rj' is reasonably good (see fig. ||) 
whereas for the a — fo no mass could be extracted 
yet. 

Gluino-glueballs. The gluino-glue correlation 
C a "(At) is evaluated by calculating the propa- 
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Table 1 



Summary of our spectrum results. New results are in bold face. 



K 


LxT 


am ++ 


am ~+ 


am. a _. n , 


amt 


arrigg 


0.1925 


12 x 24 


0.53(10) 


0.80(18) 


0.48(5) 


1.00(13) 


0.883(16) 


0.194 


12 x 24 


0.40(11) 


1.10(28) 






0.816(18) 


0.1955 


12 x 24 










0.751(21) 


0.194 


16 x 32 






0.49(6) 
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Figure 1. m e // of the ++ -glueball 
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Figure 2. m e // of the a — r/-gluinoball 



gator A"^ for a random source y. APE as well 
as Jacobi smearing for the gluino have been used 
with parameters (Nape = 9, (.ape = 0.5) and 
(iV Jac = 18, e. /ac = 0.2). 

The pseudo chiral limit. Within the OZI- 
approximation the connected term in the a — if 
correlator is expected to give rise to a massless 
mode in the limit nig — > 0, the "adjoint pion" 
a — 7r Extrapolating values for m a ^^ at var- 
ious k's, we find k c ~ 0.1970(4), which is consis- 
tent with the determination in || . 

5. Conclusions and outlook 

We still have work to do to fill out the blank 
spots in table [jj We plan to use better numer- 
ical technology such as variational smearing for 
glueballs, and stochastic estimators and spectral 
decomposition of the fcrmion propagator for the 
disconnected parts of the gluino-balls correla- 
tions. It would also be interesting to investigate 



the effects of the volume on the spectrum. For 
this better statistics are needed. 

The computations were carried out on the Cray 
T3E at NIC, Jiilich and the Sun Fire SMP- 
Cluster at RWTH Aachen, Germany. 
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